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Abstract. We perform spectroscopy on the hyperfine splitting of ^^Rb atoms trapped 
in far-off-resonance optical traps. The existence of a spatially dependent shift in the 
energy levels is shown to induce an inherent dephasing effect, which causes a broadening 
of the spectroscopic line and hence an inhomogeneous loss of atomic coherence at a 
much faster rate than the homogeneous one caused by spontaneous photon scattering. 
We present here a number of approaches for reducing this inhomogeneous broadening, 
based on trap geometry, additional laser fields, and novel microwave pulse sequences. 
We then show how hyperfine spectroscopy can be used to study quantum dynamics of 
optically trapped atoms. 



1. Introduction 

The development of laser cooling and trapping techniques 12 E] has been motivated, 
to a great extent, by their potential use for precision measurements, of which precision 
frequency measurements are of outmost importance. Cold atoms are attractive for these 
measurements mainly because of the suppression of the Doppler broadening and because 
they allow for a longer interrogation time and hence a narrower linewidth. 

Modern atomic frequency standards are based on atomic clocks, which measure 
the frequency associated with the transition between two ground state hyperfine levels. 
The performance of such a clock can be improved by lengthening the time over which 
such a transition is measured. During the 1950 's various possibilities were suggested to 
lengthen this interrogation time: Ramsey 0] proposed confining the atoms in a storage 
ring with an inhomogeneous magnetic field, or a storage cell with properly coated walls. 
Zacharias jS] proposed to build an atomic "fountain" in which a thermal beam of atoms 
is launched upward, so slow atoms from the tail of the velocity distribution will turn 
around and fall back under the force of gravity. Using the Ramsey method of separated 
fields, the atoms could be interrogated once on their way up and then again on their way 
down, resulting in a long total interrogation time. Back then, this idea failed because of 
collisions of these slow atoms with the faster ones. Using laser-cooled atoms, Zacharias' 
idea was revived in 1989 by Kasevich et. al. [H]. They cooled sodium atoms using 
a magneto-optical trap [Zj and launched them on ballistic trajectories by applying a 
short pulse of resonant light. Two short ^ microwave pulses were applied as the atoms 
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turned around inside a waveguide and allowed the ground-state hyperfine splitting to 
be measured with a linewidth of 2 Hz. 

To further increase the measurement time beyond the practical limit imposed by 
the height of a fountain, one can use trapped atoms, and in particular optical dipole 
traps, based on the dipole potential created by the interaction of a neutral atom and 
a laser beam IHlEl- Unfortunately, the interaction between the trapping light and the 
trapped atoms is not negligible. First, even for a far-off-resonance trap, where many of 
the undesired effects of the atom-light interaction are suppressed, residual spontaneous 
scattering of photons from the trap laser destroys the ground state spin polarization 
or spin coherence ^01- Next, and more importantly, atom trapping relies on inducing 
a spatial inhomogeneous shift of the atomic energy levels, and the existence of such 
a shift makes precision spectroscopy of trapped atoms a difficult task. In the case of 
hyperfine spectroscopy in optical traps, the inhomogeneity is approximately equal for 
both states, since the trapping light interacts with both levels with the same strength 
(the dipole matrix elements are identical). However, the difference in detuning between 
them results in a differential inhomogeneity, which induces a dephasing effect and, 
when ensemble- averaged, causes a broadening of the spectroscopic line and hence a loss 
of atomic coherence at a much faster rate than the spontaneous photon scattering rate 
[TT] . These effects can be reduced by increasing the trap detuning [TTl [T21 CHI E] and 
also by using blue-detuned optical traps, in which the atoms are confined mainly in the 
dark JT] E] . However, the residual frequency shifts are still the main factor limiting 
the coherence of the trapped atoms, and hence the use of dipole traps for precision 
spectroscopy. 

A number of approaches for reducing the inhomogeneous broadening of the ground 
state hyperfine splitting of optically trapped atoms are presented in this tutorial. 
Following a brief review on light shifts and optical dipole traps in section |21 and a 
description of the basic experimental setup in section IHl we present in section IH a 
simple, classical model for microwave spectroscopy in an optical trap, and use it to 
analyze the atomic coherence times for different trap geometries in section El Based 
on this model, we present in section IHl a novel method for reducing the inhomogeneous 
frequency broadening by the addition of a weak light field, spatially mode-matched with 
the trapping field and whose frequency is tuned in-between the two hyperfine levels. In 
section [7| we show another way of achieving long coherence times in an optical trap: a 
narrow energy distribution is achieved using a microwave "pre-selection" pulse which 
selects a subset of the atomic ensemble. Our technique allows the selection of a narrow 
energy band around any central energy enabling us to maximize the number of selected 
atoms (for a given energy width) by choosing the energy with the highest density of 
populated states. 

A full, quantum-mechanical model is developed in section |H1 The limit of short 
and strong pulses is further elaborated in section El and used to predict the Ramsey 
decoherence time of a thermal ensemble as a consequence of dephasing. In section 
ITUl it is shown that, for certain trap parameters, the dephasing can be reversed by 
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stimulating a "coherence echo". The failure of the echo for other trap parameters is 
due to dynamics in the trap, and thereby "echo spectroscopy" can also be used to study 
quantum dynamics in the trap even when more than 10^ states are thermally populated, 
and to study the crossover from quantum dynamics to classical dynamics. We then show 
(in section IT^ that the decay in the hyperfine coherence due to interactions with the 
environment, is only partly suppressed by echo spectroscopy, primarily due to dynamical 
(time-dependent) dephasing mechanisms. An improved pulse sequence is demonstrated, 
containing additional vr pulses for which the decay of coherence is reduced by a factor 
2.5 beyond the reduction offered by the simple echo scheme. 

Finally, in section El we show that our echo spectroscopy methods enable us to 
study quantum dynamics of trapped atoms also for chaotic and mixed dynamics where, 
surprisingly, partial revivals of atomic coherence occur in the perturbative (quantum) 
regime, and disappear in the non-perturbative (semi-classical) regime, indicating a clear 
quantum to classical transition as a function of the perturbation strength. 

2. Light Shifts and Dipole Traps 

Throughout the work presented in this tutorial we use far-off-resonance optical dipole 
traps, which are traps based on the dipole potential created by the interaction of a 
neutral atom and a laser beam [H]. A few concepts on which the rest of the work 
presented here is based, are reviewed. In sections l2 . II and (2 . 21 a brief review of the origin 
of the dipole potential, and its connection with the ac Stark shift of the ground-state is 
presented using an idealized two- level atom. In section 12.31 a short review on the use of 
the dipole potential to trap neutral atoms is given§. In sections EH and EHl a description 
of the main mechanisms for atomic decoherence is given, this time for multi-level atoms. 

2.1. The Optical Dipole Potential 

The interaction of a neutral atom and a nearly-resonant electro-magnetic field is 
governed by the dipole interaction, and is usually separated into two terms which 
correspond to a reactive force and a dissipative force. When an atom is exposed to 
light, the electric field component induces a dipole moment in the atom, oscillating at 
the driving-light frequency. 

If we assume for example a monochromatic light field with an electric field given 
by E{f,t) = eE{f)e~'^'^^^ + c.c, where e is the unit polarization vector and ujl the 
frequency, then the dipole moment is d{f,t) = ed{f)e~^^^^ + c.c. The amplitude of the 
dipole moment is related to the amplitude of the field by d = aE, where a is the atomic 
complex polarizability, which is a function of the driving frequency. The interaction of 
the induced dipole moment with the driving field gives rise to the potential [Oj 

Udiv = -\ {dE) = --i-Re(a)/(f), (1) 

§ Laser-cooling of the atoms is a prerequisite for dipole trapping, and is done with techniques that have 
became standard in many research laboratories, and will not be reviewed here (See Refs. |16lll7l[TH] l. 
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where I = 2eQc\Ef is the hght intensity. The reactive "dipole" force is a conservative 
one, and is equal to the gradient of the above potential The dissipative force is related 
to the power the oscillating dipole absorbs from the field, which is given by: 

P,,ss = (dE) = ^Im(a)/(r-). (2) 
\ / eoc 

In a quantum picture, the dipole force results from absorption of a photon from the 
field followed by stimulated emission of a photon into a different mode of the laser field. 
The momentum transfer is the vector difference between the momenta of the absorbed 
and emitted photons. The dissipative component has its origin in cycles of absorption 
of photons, followed by spontaneous emission, in a random direction. Using equation |21 
we can write an equation for the rate of spontaneous photon scattering: 

The atomic polarizability can be calculated by using the solutions of the optical 
Bloch equations, while translational degrees of freedom are taken into account |19j . 
For a two- level atom with resonance frequency cuq, and using the "rotating wave 
approximation" , the average dipole potential can be written as 

The resulting expression for the average scattering rate is 

^ i (l + ///o+V/7)0 ' 
where 7 is the natural linewidth of the atomic transition, 5 = ujl — ujq\s the detuning of 
the laser from the atomic resonance, Jq is the saturation intensity of the transition, given 
by Jo = 27r^%c/3Ag, and Aq = 2t:c/ujq. For a large detuning from resonance (5 ^ 7) 
equations m and El can be approximated as 



2.2. Ground State Light Shifts 

It is useful to look at the atom-light interaction in the "dressed state" model ^H]- The 
combined Hamiltonian for the atom and the laser field is given by 7i = TIa + Ti.L + Ti-AL, 
where TYa^^l, and TYal are the atom, laser and interaction parts of the Hamiltonian, 
respectively. Consider a two- level atom with ground state \g) and excited state |e) with 
energy hwo, in the presence of a laser with frequency ul- The atom can be described 
by the Hamiltonian 

Ha = hwo \e) (e| , (8) 
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Figure 1. Dressed state picture of a two- level atom coupled by a light field with 
frequency below the atomic resonance, (a) The energy levels form a ladder of (almost 
degenerate) manifolds separated by the photon energy hiUL- The atom- field interaction 
separates the levels within each manifold by h£l' , where fi' is the effective Rabi 
frequency. The eigenstates of the coupled system are a mixture of the eigenstates 
of the uncoupled system, (b) If ^ \5\, then the dressed eigenstates can be identified 
with the original ground and excited states, with a small portion of the other state. The 
energy shift is Ml"^ /AS. (c) A spatially inhomogeneous light field produces a ground- 
state potential well, in which the atom (shown schematically for a Gaussian trap) can 
be trapped. 



and the laser field by 



nL = hwL(a^a+]-] , (9) 



2^ 

where and a are the photon creation and annihilation operators, respectively. If the 
coupling between the atom and the field is ignored, then the eigenstates of TIa + 'Hl 
are characterized by the atomic internal state (either \g) or |e)) and by the number of 
photons in the electromagnetic field, n. The energy levels form a ladder of manifolds, 
separated by huJi, each containing two states of the form \g,n) and \e,n — l) (see 
figure ^). If the light is in resonance with the atomic transition, these two states 
are degenerate. Otherwise, they are separated by h\5\. 

A general result from second-order time-independent perturbation theory is that 
an interaction Hamiltonian Ti' will lead to a shift in the energy of the (unperturbed) 
i-th state given by |9J 

^^.-E (10) 



The interaction term in our case is TYal = — d ■ i?, where d = — er is the electric dipole 
operator. This interaction term couples states within the same manifold with a matrix 
element given by 

{e,n-l\HAL\g.n) = —, (11) 

where Vt = dE/h is the Rabi frequency. The eigenfunctions of the coupled system, 
denoted "dressed states", are a mixture of the uncoupled system eigenstates and are 
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given by [TH] 

|l(n)) = cos 6' |e,n — 1) — sin 6' \g,n) (12) 
\2{n)) = sin 6* I e, n — 1) + cos^^ \g,n) , (13) 

where tan2^^ = —Q/6. The new (dressed) levels are separated by an energy HQ', 
where Q' = y/Q'^ + 6"^ is called the "generalized Rabi frequency" . For a large detuning, 
\6\ ^ Q, the ground state is shifted by (see figure [Dd). The dressed states 

|2(n)),|2(?2 + 1)), ... can be identified with the original ground state, with a small 
mixture of the excited state. This light-induced shift in the atomic energy levels is 
usually denoted "light shift" or "ac Stark shift". Using fi^ = (7^/2) (I/Iq), where Iq is 
the saturation intensity of the transition, we conclude that the light-shift of the ground 
state energy, AEg, is given by 

As seen, in the perturbative limit the light-shift of the ground state is equal to the dipole 
potential in equation IHl The reason is that, since in this regime the atom is mostly in 
the ground state, the light shifted ground state can be identified as the potential for 
the motion of the atoms (see figure^). In the case of multi- level atoms, equation IT^ 
should be modified to include the electric dipole interaction between the ground-state 
and all the excited states, with their respective detuning and transition strength (see 
section 12. 5p . 

2.3. Optical Dipole Traps 

The dissipative part of the atom-light interaction is used for laser cooling of atoms 
[T^ . a pre-requisite for optical trapping. However, spontaneous photon scattering is in 
general detrimental to trapped atoms, mainly because it can induce heating and loss. 
Comparing equations IHl and [7| yields. 



Udip S 



(15) 



his 7 

which results in the well known fact that a trap with an arbitrarily small scattering rate 
can be formed by increasing the detuning while maintaining the ratio 1/6. 

Equation IHl indicates that if the laser frequency is smaller than the resonance 
frequency, i.e. 6 < ("red-detuning") the dipole potential is negative and the atoms are 
attracted by the light field. The minima of the potential is found then at the position 
of maximum intensity. In the case 6 > { "blue detuning" ) the minima of the potential 
is located at the minima of the light intensity. 

Trapping atoms with optical dipole potentials was first proposed by Letokhov [20] 
and Ashkin fIT\ . Chu and coworkers jH], were the first to realize such a trap, trapping 
about 500 atoms for several seconds using a tightly focused red-detuned beam. Later, 
a far-off- resonant trap for rubidium atoms was demonstrated ^3], with a detuning of 
up to 65 nm, i.e. 6 > 5 x 10^7. In this case, the potential is nearly conservative and 
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spontaneous scattering of photons is greatly reduced. A comprehensive review of the 
different schemes and apphcations of such optical dipole traps is presented in reference 
[Oj. In the limiting case where the frequency of the trapping light is much smaller then 
the atomic resonance, trapping is still possible, practically with no photon scattering 
[22] • Such a quasi-electrostatic trap, formed by two crossed CO2 laser beams, was used 
to create a Bose-Einstein condensate without the use of magnetic traps [22] • 

Apart from using far-off-resonance lasers, the interaction between the light field and 
the atoms can be reduced by the use of blue-detuned traps, in which atoms are confined 
mostly in the dark J3] . Experimentally, blue-detuned traps, which require surrounding 
a dark region of space with a repulsive dipole potential, are harder to realize than 
red-detuned ones, where already a single focused beam constitutes a trap [H]. Several 
configurations for far-detuned dark optical traps were demonstrated, in which gravity 
provided the confinement in one direction: "Light sheets" traps were generated by 
elliptically focusing two laser beams and overlapping the two propagating light sheets to 
form a "V" cross-section that supports against gravity, while the confinement in the laser 
propagation direction is provided by the beam divergence ^1]. In a later work, and in 
order to achieve larger trapping volume, a different trap was constructed with four light 
sheets producing an inverted pyramid 24^ . A single beam trap was demonstrated using 
two axicons and a spherical lens to generate a conical hollow beam propagating upwards 
Such gravito-optical traps are limited to weak confinement. Traps in which light 
provided the confinement in all directions were developed with hollow beams. Laguerre- 
Gaussian "doughnut" modes were used, together with additional plug-in beams, to form 
such a trap ^B] . 

Several dark traps based on a single laser beam were demonstrated, providing 
grater experimental simplicity and enabling dynamical changes of the trap geometry 
and strength. As opposed to the 2D case, were any desired light distribution can be 
generated using diffractive or refractive optical elements, there is no simple procedure 
to design an arbitrary 3D light distribution [2Z1 12H] • Nevertheless, through the use of 
refractive and holographic optical elements, it is possible to produce light distributions 
which are suitable for trapping atoms in the dark using a single laser beam. Such light 
distributions, which comprise of a dark volume completely surrounded by light, were 
reahzed using either combinations of axicons and spherical lenses, diffractive optical 
elements, or rapidly scanning laser beams. 

In the first scheme realized, the trapping beam was produced by passing a Gaussian 
beam through a phase plate of appropriate size, which shifted the optical phase at the 
center of the beam by vr radians. Interference leads to a dark volume at the focus of 
the lens, surrounded by light in all directions . An additional method, with a much 
larger volume and more symmetric shape, was realized by simultaneously focusing two 
diffraction orders of a properly designed binary phase element, consisting of concentric 
phase rings with a vr-phase difference between subsequent rings [20] • An improved trap 
configuration was demonstrated by adding an axicon telescope before the phase element 
of the above setup |31J. This configuration maximizes the trap depth for a given laser 
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power and trap dimensions, and greatly reduces the light induced perturbations to the 
trapped atoms. A related scheme was demonstrated in references (32] and ^3] ■ Finally, 
a tightly- focused rapidly-rotating laser beam was used to create a trap |31j • If the scan 
frequency is high enough, the optical dipole potential can be approximated as a time 
averaged quasi-static potential. For a blue-detuned laser beam, and a radius of rotation 
larger than the waist of the focussed beam, a dark volume suitable for 3D trapping is 
obtained. 



2.4- Photon Scattering and Coherence Relaxation in Multi-Level Atoms 

The dipole potential can be viewed as originating from cycles of absorption of a photon 
and stimulated emission into a different laser mode. This process is unfortunately 
accompanied by spontaneous scattering, in which the absorption is followed by 
spontaneous emission in a random direction. Spontaneous scattering is one of the main 
limiting mechanism for the atomic hyperfine coherence||. 

Spontaneous scattering is a two-photon process, in which an atom initially at a 
state |F, ttt-f) absorbs a photon from the trap laser and moves to an intermediate state 
\F', m'p) of some excited level. The atom then decays back to the ground state, to a final 
state \F'\m"p). If F" = F and m'p = mp, the process is called "Rayleigh scattering". 
Otherwise, it is denoted "Raman scattering" . Clearly, Raman events destroy the atomic 
state coherence and hence are detrimental to the spectroscopy. As will be discussed 
below in section [T21 also Rayleigh scattering events can be destructive for spectroscopy 
of a trapped ensemble. 

The probability amplitude for scattering between \F,mF) and \F",mp), via an 
intermediate state \F',m'p) is proportional to {F",m'p F',m'p) {F',m'p F,mp), 
where fik if^j) are the spherical components of the dipole moment operator {k,j = —1, 
or 1 depending on the polarization of the absorbed and emitted photons). If the detuning 
from resonance is large enough such that no specific intermediate state is resolved, 
then to calculate the rate of transitions 'yF,mF^F",m'^ between |F, mj?) and \F",m'p) 
the amplitudes for all possible paths must be summed. For rubidium atoms, these 
are the two excited states and 5^P3/2, with all their different hyperfine levels and 

corresponding Zeeman sub-levels. We further assume that the laser polarization is linear, 
hence conservation of angular momentum dictates mp = m'p. The total transition rate 
is HO] 



F 



Sttc'^ujII 
2/i/i4 



(1) (2) 
^F,mF^F"m'^ ^F,mF^F"m'^ 



Si 62 



(16) 



II A word of caution is in order about our use of the terms "decoherence" and "dephasing" . Decoherence 
results from the dissipative interaction of a single quantum superposition state with the environment 
(e.g. in the case of trapped atoms, the coupling to the electro-magnetic vacuum which leads to 
spontaneous scattering of photons). In addition, the response of a macroscopic ensemble of quantum 
systems (each prepared in a superposition state) decays due to the dephasing between the microscopic 
systems resulting from local variations in their evolution, and hence the ensemble coherence is lost. 
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Figure 2. Calculated total scattering rate (dashed line) and i^-changing Raman 
scattering rate (full line) for ^^Rb atoms and linear polarization. The former affects the 
coherence time of trapped atoms, but the later is experimentally simpler to measure. 



where 

o^fIip^f" m''=-^ ^"f I/^?! F'^ ^'f) {F\ ^'f I/^oI F, THf) ■ (17) 

Here ujj and 7j are the transition frequency and excited state hfetime of the Dj hne 
(J = 1, 2), and /i = (33 |/i-i| 44) is the amphtude for the strongest transition. 

The first step in spectroscopy of trapped ®^Rb atoms, is to transfer the entire 
population to F = 2. There, it is distributed uniformly between niF levels. We then 
apply a certain microwave pulse sequence (e.g. a Rabi pulse) which acts only on the 
transition \F = 2,mF = 0) —>■ \F = 3,mF = 0), since a bias magnetic field is applied 
to shift all other transitions from resonance. Finally, we detect the total population 
of F = 3. We would like to measure, independently, the rate of incoherent scattering 
events that also affect the population of F = 3, in order to subtract them from our 
coherently driven population changes. Experimentally, we can measure the amount of 
F-changing transitions, but not the rate of Rayleigh scattering events, or Raman events 
that do not change F ^0] • Equation ^1 together with the appropriate summation, is 
used to calculate the relative amount of different scattering processes. For example, 
the rate '-fF=2^F=3 for a transition between hyperfine levels 2 and 3 is calculated by 
averaging '^F^mp^F" ,m'^ over initial mp levels of F = 2, and summing over final m"p 
levels of F = 3. The total scattering rate is given by summing over mp levels of F = 2 
and m"p levels of F = 3. Figure |21 shows a calculation of the relative amount of these 
two processes, as a function of the laser wavelength. Using the results of figure El we 
calculate the total amount of scattering events out of our measurement. 

Note, that for a detuning larger than the fine-structure sphtting of the excited state 
(15 nm in rubidium), a destructive interference exists between the transition amplitudes 
for Raman scattering, summed over the intermediate excited states p^. In this case. 
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most spontaneous scattering events leave the internal state of the atom unchanged. 



2.5. Light Shifts and Dephasing in Multi-Level Atoms 

For a real multi-level atom, the light-shifts described in section 12.21 depend, in general, 
on the particular substate of the atom. The equivalent of equation El in the case of a 
multi-level atom can be written using the dipole matrix elements fiij = {ci \fi\ gi) between 
a specific ground state \gi) and a specific excited state \ej). Using fiij = Cij we can 
write the ac Stark shift for a ground state \gi) as [H]: 

The summation takes into account the contributions of the different coupled excited 
levels \ej), each with its respective transition coefficient Qj, and detuning 6ij = u — Uij. 

If the detuning of the light is large as compared to the excited state hyperfine 
splitting cij^p, then certain sum rules exist for the transition coefficients [SH] and equation 
^Ifor the shift of a ground state with total angular momentum F, simplifies to: 

TTc'^'yl 



^ + Y~) - gprnpy/l 

h,F (>2,F / 



e2 



(19) 



where e is the light polarization ellipticity and 6i^f [^2,f] is the detuning of the laser 
from the Di [D2] line. For linearly polarized light (e = 1) 

AE.(r) = !j^IM. (20) 



where 



1 _ 1 

si's 



A + A ' A (21) 

is the "weighted detuning" (below, we drop the "*" and call this the detuning). Note, 
that the interaction matrix elements are identical for atoms in \F = 2) and \F = 3). 
However, the light shift is inversely proportional to the trap laser detuning, which differs 
by t^HP, and therefore the potential from a far-detuned laser is slightly different for the 
two hyperfine levels. Loosely speaking, this difference in potential means that different 
hyperfine states "feel" different trap shapes. For microwave spectroscopy of a thermal 
ensemble this results in a rapid decay of the ensemble-averaged spectroscopic signal as 
a consequence of dephasing, much faster than the decay due to spontaneous scattering 
of photons. 



3. Experimental Setup 

All the experiments described in this tutorial consist of three stages. During the first 
stage standard laser-cooling techniques are used to cool an ensemble of ^^Rb atoms 
and load a portion of them into a far-off-resonance trap. During the second stage, all 
the nearly-resonant laser beams are shut-off leaving only the far-off-resonance trapping 
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beam on, and microwave pulses are applied. In last stage short pulses of on-resonance 
beams are used for the diagnostics. In this section, a detailed description of the above 
steps is provided. 

3.1. Laser Cooling and Trapping 

The heart of our setup is a vacuum chamber connected to a small reservoir of rubidium 
atoms through a valve. The chamber has six big windows used for the laser beams, 
and two additional small windows used for imaging the atomic cloud with an intensified 
CCD camera and measuring the fluorescence signal with a photomultiplier tube (see 
figure EI) . 

The first step in all the experiments described here is a magneto-optical trap (MOT) 
(Tj. The MOT laser beams consist of three orthogonal pairs of counter-propagating 
beams, detuned approximately —17 MHz (= —2.87) from the 155*1/2,-^ = 3) — >• 
I5P3/2, -F' = 4) line, and a "repumping" beam in resonance with 155*1/2,-^ = 2) — »• 
I5P3/2, -F' = 2). A pair of water-cooled copper coils in the anti-Helmholtz configuration 
provide the magnetic field gradient for the MOT (see figure 121). In addition, three 
orthogonal Helmholtz coils are used to compensate for constant magnetic fields. 

The MOT loading time is ~ 700 ms. After that a ~ 50 ms temporal dark MOT 
stage [12] is applied in order to increase the spatial density of the atoms. The intensity of 
the MOT beams is reduced to ~ 1.5 mWcm~^, their detuning increased to ~ —30 MHz, 
the magnetic field gradient increased to 8 Gcm~^, and the repump intensity reduced by 
a factor of 40. Final optimization of these parameters is performed by optimizing the 
number of atoms loaded into the dark optical trap. A peak density of ~ 1 x 10^^ cm~^ 
is achieved for optimized parameters, but typically we work at lower densities, for which 
collisions are negligible for the timescales of our experiments. The laser cooling stage 
ends with a polarization-gradient cooling stage [HH]; in which the magnetic field is set 
to zero, the MOT beams intensity are further decreased and their detuning increased. 
Usually, 3 ms of this cooling stage result in a temperature of ~ 10 — 20 /iK. 

The dipole trap in our experiments consists of a linearly polarized horizontal 
Gaussian laser beam with wavelength in the range A = 785 — 810 nm. The trap beam is 
coupled into a single mode fiber, and the output of the fiber (with power in the range 
P = 50 — 400 mW) is focused to a ~ 50 fim spot at the center of the vacuum chamber. 
An active servo circuit and an acousto-optic modulator ensure a 1% stability in the beam 
power, and hence in the trap depth. With typical values of P = 50 mW and A = 800 nm 
we achieve a trap with a depth of ~ 35 /iK, and oscillations frequencies Ur = 2.3 KHz 
and uJz = 8.4 Hz in the radial and axial dimension, respectively. The clear separation of 
time-scales between the fast transverse oscillations and the very slow longitudinal ones 
is used in the analysis of some of our experiments to neglect the longitudinal motion, 
which is essentially frozen during the duration of the experiment, and hence to treat the 
system as a two-dimensional one. 

The trap beam overlaps the center of the atomic cloud during the cooling stages 
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Figure 3. Basic components of the experimental setup. The MOT beams consists 
of three orthogonal pairs of counter-propagating beams, and a repumping beam. Two 
copper coils in the anti-Helmholtz configuration provide the magnetic field gradient, 
and three orthogonal Helmholtz coils are used to compensate for constant magnetic 
fields. The atoms are imaged with an intensified CCD camera, and their fluorescence 
measured with a photomultiplier tube. The dipole trap beam enters the chamber, 
co-propagating with the detection beam. The microwave pulses are created by a 
synthesizer, amplified and radiated into the chamber using an antenna. A bias 
magnetic field is applied parallel to the trap's polarization axis and to the microwave 
magnetic field direction, in order to Zeeman shift the magnetic sensitive mp ^ levels 
out of resonance with the microwave pulse. 



and, after turning off all the lasers (with the exception of the trap) we end up with 
~ 10^ atoms loaded in the trap (depending on the specific power and detuning) with a 
temperature of ~ 20 fiK. At the end of this stage, the atoms are prepared in the F = 2 
ground state by turning on the MOT beams, without a repump beam j^HJi for 1 nis. 

3.2. Microwave Spectroscopy 

We are interested in spectroscopy of the "clock" transition, i.e. the cUhf = 27r x 3.032 x 
lO^sec"^ transition between the two magnetic insensitive hyperfine Zeeman substates 
of the ground state, \F = 2,mF = 0) and \F = 3, mi? = 0), which we denote |J.) and 
It). We drive this transition with a nearly-resonant microwave field. The microwave 
pulses are created by a synthesizer (Anritsu, 69317B), whose clock is connected to a 
high stability and low phase noise 10 MHz oscillator. The pulses are then amplified 
with a 30 dB amplifier (Mini-Circuits, ZVE-86) and irradiated into the chamber using a 
Log-periodic antenna. The strongest microwave fields produced in our setup correspond 
to a Rabi-frequency of 5 kHz for free (untrapped) atoms. 

A bias magnetic field is applied parallel to the trap's polarization axis and to the 
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microwave magnetic field direction, in order to Zeeman shift tlie magnetic sensitive 
nip ^ levels out of resonance with the microwave pulse. In most of the experiments 
(See sections El and E]]) 5 its value is ~ 40 — 80 mG. For the experiment in section IT^ it 
is ~ 240 mG, but then a MOSFET switch is used to turn it on after the cooling stage, 
which requires a nearly zero field. Special care is taken to align the directions of the 
microwave magnetic field and the bias magnetic field to be parallel to the trap's laser 
magnetic field (i.e. the trap's polarization is perpendicular to the bias field), to enable 
a common well-defined quantization axis for all the fields interacting with the atoms. 

3.3. Diagnostics 

Following the microwave pulses, (the population in F = 3) is measured by detecting 
the fluorescence during a short pulse of a laser beam resonant with the cycling transition 
155*1/2, -F = 3) ^ |5P3/2,-F = 4). The population of F = 2 is then measured for the 
same experimental run by turning on the repumping beam (which is resonant with 
|5S'i/2, -F = 2) — » |5P3/2,F = 3)) and applying an additional detection pulse. This 
normalized detection scheme is insensitive to shot-to-shot fluctuations in atom number 
as well as fluctuations of the detection laser frequency and intensity j37j. 

In addition to the population of |t) by the microwave field, ||) and all sublevels of 
F = 3 are in general populated by spontaneous scattering of photons from the trapping 
laser, which always accompanies the dipole potential and tends to destroy the atomic 
coherence (see sectional). A simple way to measure the amount of photon scattering 
is by measuring the spin relaxation caused by spontaneous Raman scattering [TUl I29j. 
This is a very useful experimental technique which enables to measure even very low 
scattering rates. For this measurement, the trapped atoms are first prepared in the 
lower hyperfine level of the ground state, F = 2. The number of atoms in F = 3 after 
a variable time t, (t), is measured by detecting the fluorescence after a short pulse 
of a resonant laser beam. The fraction of Raman scattering events out of the total 
photon scattering events can be calculated using the known matrix elements for the 
transition, and therefore the total amount of photon scattering can be inferred from the 
spin relaxation measurement (see section l2^ . 

We calculate the microwave driven part of the population by subtracting from the 
measured signal the above contribution due to F-changing Raman transitions induced 
by the trap laser and normalize to the signal after a short vr-pulse, which transfers the 
whole population of ||) to ||). The corrected and normalized signal is denoted Pp 

4. Microwave Spectroscopy in a Dipole Trap: A Classical Model 

We drive the magnetic-insensitive transition with a nearly-resonant microwave field, and 
apply a static magnetic field in order to shift all other (magnetic sensitive) levels out 
of resonance. In such case, and neglecting for the moment the effects of the trap on 
the atoms, we can consider a two- level system, separated in energy by Fhf = /ao^hp- 
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Figure 4. Rabi spectrum of free atoms with a 25 ms microwave pulse, showing the 
characteristic "sinc^" hneshape corresponding to a rectangular pulse. Note that since 
the population of the four |F = 2, m_F ^ 0) states is included in iVf / {Ni +Ni), a value 
of 0.2 represents the maximal possible signal (a tt pulse) for the \F = 2,mF = 0) ^ 
l-F = 3, mi? = 0) transition. 



When a free (untrapped) atom, initially in state |J,), is placed in a microwave field with 
frequency cJmw the probability of finding it in state |t) as a function of time is: 

PT = ^sin^(%^t), (22) 

where f^Mw is the Rabi frequency of the microwave field, ^'^^ = a/^^mw + the 
generalized Rabi frequency, A = cJmw — ^^hf is the microwave detuning, and we have 
used the rotating wave approximation. As a function of the interaction time, an atom 
oscillates between the two states at the effective Rabi frequency, where the contrast of 
the oscillation is determined by the Rabi frequency and the detuning. This sinusoidal 
population transfer is referred to as Rabi flopping. If the interaction time is fixed and the 
power is varied, reaches a maximum for flmwi = vr. For such a pulse, the transition 
probability clS db function of frequency has a maximum for cUmw — '^hfi 

and (FWHM) 

width of Auj/2'K ^ 0.8/t. A measurement of the resonance frequency using a single 
pulse of constant duration is called Rabi spectroscopy. Figure IH shows a Rabi spectrum 
of cold atoms freely falling after the shut off of the cooling and trapping beams. The 
pulse duration is limited by the falling time of the atoms from the detection volume. For 
this short pulse we observe a Fourier-limited linewidth, indicating that no decoherence 
or dephasing mechanisms exist for for free-falling atoms in this time-scales. 

An alternative and widely used method is Ramsey spectroscopy, in which two 
consecutive pulses of duration t, separated by a "dark" period of duration T, are used. 
The first pulse (called a | pulse) puts the atoms in a quantum superposition of both 
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states. During free evolution, due to the detuning of the driving field from resonance, 
these two states develop a relative phase in the rotating frame of the field. The second | 
pulse causes the two states to interfere, thus translating phase differences into population 
differences. If the field frequency is scanned, the phase difference acquired in the dark 
period, and thus also the population of the levels after the second pulse, oscillate. 
Following a two-pulse Ramsey sequence, the probability of being in the excited state as 
a function of time for |A| ^ ^mw and t <^ T is simply 

A = ^ + ^ cos (AT) . (23) 

When performing Rabi or Ramsey Spectroscopy on optically trapped atoms, the 
above treatment no longer holds. The dipole potential is inversely proportional to the 
detuning of the trapping laser beam from resonance. Since there is a slight difference in 
potential for atoms in different hyperfine states, the external (center of mass) potential 
depends on the internal (spin) state, hence the internal and external degrees of freedom 
cannot be separated. 

We first consider a simple model of Rabi spectroscopy, in which we assume the 
atoms to be in thermal equilibrium in the potential of the trap and frozen in position 
during the microwave pulse. In order for this model to be valid, the microwave pulse 
duration t must be much shorter than the typical oscilation time in the trap, i.e. 

t « (24) 

Let us look now at the ground states of atoms trapped in a dipole trap. From 
equation 1201 it is seen that in the presence of the trap light (and assuming for simplicity 
that 62,F ^ Ai?, so the contribution from the Di line is neglected) the energy splitting 
is modified by 



huj^p (r) - huJuF = — 3 — TT 
UJq 



9 -K-j, (25) 

where uJ^ (r) is the spatially dependent transition frequency in the presence of the 
light and 5 = {52,f=2 + 52,f=3)/2 measures the laser detuning from the center of the 
ground state hyperfine splitting. Equation |2S1 indicates that for \6\ > i.e. for 62,f=2 
and S2,F=3 both positive (or both negative), the ground state energy splitting is always 
reduced by the presence of a light field. When the detuning is "between" the hyperfine 
levels, \S\ < the energy splitting is enlarged. 

For a trap with detuning 5 ^ equation |2S1 yields 

Uhf (r) - a;HP ~ j — (26) 

where U{r) is the spatially dependent dipole potential that forms the trap. The fact 
that the relative ac Stark shift is o^hf/ S smaller than the dipole potential, is the main 
motivation for using far-off-resonance traps for precision spectroscopy. For example in 
reference jTT] the relative ac Stark shifts were only u^p/'S ~ 2 ■ 10~^ times the dipole 
potential. 
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For a trapped atomic ensemble, averaging over the different values of U{t) causes 
a shift in the ensemble averaged transition frequency , {lJ^ — uj^p). This shift is also a 
function of the trap's geometry, which determines the degree of exposure of the atoms 
to the trapping light. Following references jH] and ^3], we introduce the parameter 
K, defined as the ratio of the ensemble-averaged potential and kinetic energies of the 
trapped atoms, k = (U) / {Ek), and refer to it as the "darkness factor" of the trap. 
Assuming a trapped atomic gas in thermal equilibrium, the ensemble- averaged kinetic 
energy is {Ek) = ^ksT, and neglecting gravity, the ensemble averaged potential energy 
is given by 



J drU (r) exp 


\ f/(r)-f/ol 




J dr exp 


\ Uir)-Uo] 






ksT 





(27) 



where Uq is the potential at the trap bottom and the integration is over the entire trap 
volume^. The ensemble frequency shift is then given by 

{uj^ - Uiip) = ^ • ^kbT ■ K. (28) 

Using equation [3 we can calculate the average spontaneous scattering rate, which 
determines the homogeneous coherence time of the atomic ensemble: 

(7s) = Ys ■ ■ (29) 

In addition, the spatial dependence of U{r) will result in an inhomogeneous broadening 
responsible for the inhomogeneous decoherence time: 



aM ^ y/(^^> - i^pf = ^^{Uirf) - (t/(r))l (30) 
For example, for a thermal ensemble in an harmonic trap, we have 

(7.) = i (IhT) , (31) 



hS \2 
and 

Using equations ED and 1221 we can see that 

a{uj^p) ujp_ 



. (33) 

{is) 7 

For most alkali-metal atoms, such as rubidium and cesium, uJuf/i ~ 10'^, indicating 
that, however small, the relative ac Stark broadening is still much larger than the 

% For red-detuned traps, Uq < 0, while Uq — for most blue-detuned traps. For example, for an 
harmonic trap, the ensemble-averaged potential energy (relative to the trap's bottom) and the averaged 
kinetic energy are equal, and therefore the darkness factor is always k = 1 for a blue-detuned trap, 
independent of laser power, detuning or atomic temperature, but k cx Uo/ksT for a red-detuned trap 
with Uo » ksT. 
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spontaneous photon scattering rate. The homogeneous and inhomogeneous coherence 
times are inversely proportional to the spontaneous scattering rate and the ac Stark 
broadening, respectively, hence equation IHSl shows that the latter is the main limiting 
factor for the atomic coherence time in the trap. 

Note, that in order to resolve the frequency distribution, the Fourier broadening 
must be smaller than the inhomogeneous broadening, i.e. t~^ <^ cr(ti;HF)- Hence, using 
also equation 1211 we conclude that the requirement for the present classical "stationary" 
model to be both valid and "interesting" is 

uJosc < o-(t^Hp)- (34) 

5. Case Study: Inhomogeneous Broadening and Trap Geometry 

As a specific example for the use of the classical model presented in the previous section 
we compare the performance of four different trap geometries. One trap is red-detuned 
(the crossed Gaussian beams trap [111), ^^^1 three are blue-detuned: the Rotating Beam 
Trap (ROBOT) j24| is based on a repulsive optical potential formed by a tightly focused 
blue-detuned laser beam which is a rapidly scanned using two perpendicular acousto- 
optic scanners. The Laguerre-Gaussian (LG) trap 1221, consists of a hollow laser beam 
(LGg) and two additional "plug" beams that confine the atoms in the propagation 
direction of the hollow beam. Finally, the "optimal trap" [3T] uses two refractive axicons 
and a binary phase element to create a spatial light distribution consisting of two hollow 
cones attached at their bases and completely surrounding a dark region. 

Figure El shows the calculated optical potential distribution of three of these traps, 
in the r—z plane. The dashed lines in figureElare equidistant contours of equal potential, 
and the solid one is the contour line corresponding to the trap depth. We assume a trap 
laser with a fixed power P = IW and a sample of ^^Rb atoms, which are laser cooled 
to a temperature of 5 fiK, and form a nearly spherical cloud with a radius of 0.5 mm, 
typical parameters for a magneto-optical trap. 

Adopting a criteria of > 90% geometrical loading efficiency from the magneto- 
optical trap, we choose a radius r = 0.5 mm for all dark traps. The beam waist is 
chosen as wq = 50 fim for the ROBOT, and wq = 10 /im for the optimal trap. The 
length of the latter is an independent parameter, chosen as L = 3 mm to optimize 
the power distribution as explained in reference [SI]. We neglect the enhanced loading 
efficiency of red-detuned traps [SHI, ^^^1 assume for the crossed trap Wq = 0.6 mm, which 
corresponds to > 90% overlap with the magneto-optical trap [T^ . 

The detuning in the comparison is chosen such that the depth of each trap is 3 
times larger than the mean kinetic energy of the atoms. Since a fixed laser power 
is assumed, less efficient traps would require a smaller detuning to provide the same 
trap depth. In table El the calculated required detuning in each of the different 

+ We choose a crossed trap, and not a simpler focused gaussian beam trap, since with a single focused 
beam a trap radius of 0.6 mm will result in an extremely large axial size of > 1 m. 
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Figure 5. Contour maps of the calculated trapping depth for three different optical 
traps. The dashed lines are equidistant contours of equal potential. The solid line 

is the contour corresponding to the trap depth. All the traps have the same radial 
dimension, (a) Crossed Gaussian beam (red-detuned) trap, (b) Rotating beam trap, 
(c) "Optimal" trap. 
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6 (nm) 




(7.) (s-^) 


cr(a;HF) (Hz) 


Red detuned trap 


-0.7 


4.9 


166.9 


1.9-103 


LG trap 


0.23 


0.6 


87.6 


5-103 


ROBOT 


0.19 


0.2 


21.2 


3.8-103 


"Optimal" trap 


4.69 


0.02 


0.09 


47.2 



Table 1. Required detuning, calculated atomic darkness factor, mean spontaneous 
photon scattering rate and inhomogeneous frequency broadening of the hyperfine 
splitting for ^^Rb atoms confined in the traps analyzed in the text. 

traps, with the parameters discussed above, is presented together with the calculated 
darkness factor n. As expected, all blue-detuned traps have a better darkness factor 
than the red-detuned trap. The optimal trap has a significantly better darkness factor 
{k = 0.02) than all other schemes, due to its very thin walls and nearly minimal surface 
to volume ratio. Next, the mean spontaneous photon scattering rate (7^), and the 
inhomogeneous hyperfine frequency broadening cr(u;HF) are calculated using equations 
EHl and EDI respectively. Here, the advantage of the optimal trap is even much larger, 
since the improved darkness factor is combined with the efficient distribution of optical 
power that enables an increased detuning for the same trap depth. For example, the 
inhomogeneous broadening in the optimal trap is only 47 Hz, while for all the others is 
of the order of a few KHz. 

Including gravity in our calculations, results in an increase of 10-60% for the 
scattering rate and inhomogeneous broadening for rubidium atoms in the above dark 
traps. For the lighter alkali-metal atoms (e.g. sodium and lithium), the inclusion of 
gravity yields an increase of no more than 5% in the scattering rate. 

Our classical model and the above comparison show that, in addition to the use 
of far-off-resonance lasers and dark optical traps, geometry can also be exploited to 
minimize the undesired exposure of the atoms to the trapping light. 

6. A Compensated Trap 

The previous section indicated that even for far detuned traps with favorable geometries 
the hyperfine coherence of optically trapped atoms is still predominately limited by the 
difference in the trap-induced ac Stark shifts between the two hyperfine levels. To cancel 
these relative ac Stark shifts, we introduce an additional laser beam, with intensity /' (r) 
and frequency between the resonant frequencies of the two ground state hyperfine levels, 
say in the middle, i.e. — C(Jhp/2 and +uj-iip/2 from the lower and higher hyperfine level 
respectively*. The total shift is obtained by adding the shifts from the trap and the 

* The detuning of the compensating laser can be chosen in the range — ^^^f^ < S < yielding 
a straightforward modification of equations K^5I36I however, choosing (5 = minimizes spontaneous 
photon scattering. 
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Figure 6. Ground level energies (a) and energy difference (b) of atoms trapped in a 
focused Gaussian beam. When exposed to the trapping laser, the two hyperfine levels 
have a different ac Stark shift (dashed line). An additional weak laser, detuned to the 
middle of the hyperfine splitting, creates an ac Stark shift (dotted line) such that the 
total amount of light shift (full line) is identical for both hyperfine levels. 



"compensating" beam, 

I(r) /' 



"-•^HP (rj — CUhP — q 

^0 



,2 



(35) 



If the compensating beam is spatially mode-matched with the trap beam, i.e. /' (r) = 
?7 X / (r), then a complete cancellation of the inhomogeneous broadening will occur for 

Figure IHl shows the shift of the hyperfine levels (a) and the hyperfine energy 
difference (b) caused by the trapping beam (dashed line) and compensating beam 
(dotted line). In the presence of both beams, the levels are shifted by the same amount 
(full line). As a specific example, with a trap detuned by 5nm we have r] ~ 3.6 x 10~^. 
Hence, with a typical trap power of 50 mW, the required compensating beam power 
is 20 nW. (A similar calculation can be made which takes into account also the 
contribution of the Di transition, and introduces only a small correction to equation 



Note, that the dipole potential created by the compensating beam is U' (r) = 
±|^^[/ (r) for atoms in the upper an lower hyperfine level, respectively, and hence is 
negligible when compared with the potential of the dipole trap. Moreover, the photon 
scattering rate from the compensating beam 7^ is given by h'j'g ~ 2'-fU' /uuf, which can 
be also written as hj^ ~ jU (r) ^ h'js- Hence, the scattering rate from the nearly 
resonant compensating beam, is of similar magnitude to that of the far-off-resonance 
trapping beam, 7^. 

We implement the proposed scheme with a red-detuned (A = 785 nm) Gaussian 
trap P], created by focusing a 50 mW laser to a waist of Wq = 50 /im (See section |2I) 
and resulting in a potential depth of Uq ~ 34 fiK and oscillation frequencies of 2300 and 
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Figure 7. Rabi spectrum of the hyperfine sphtting of ^^Rb, with a 3 ms pulse, (a) 
Spectrum of free falling atoms, (b) Spectrum of trapped atoms, showing a shift in 
the line center and a broadening, (c) Spectrum of trapped atoms, with an additional 
compensating beam. The addition of the weak compensating beam, nearly cancels the 
shift and broadening of the spectrum. 



8 Hz in the radial and axial directions, respectively. An additional laser, with frequency 
locked close to the middle of the ground state hyperfine splitting, is combined with 
the trap laser jSHj- To achieve optimal spatial mode- match, both lasers are coupled 
into a polarization-preserving single-mode optical fiber, and the fiber's output is passed 
through a polarizer and focused into the vacuum chamber. Two servo loops are used to 
control and stabilize the power of the lasers: The first one ensures a 1% stability of the 
trap laser. More importantly, for complete compensation of the relative ac Stark shifts, 
a second servo loop ensures a 0.1% stability of the power ratio rj throughout the entire 
duration of the experiment. Since typically rj ~ 10~^ in our experiment, the beams 
are separated by two gratings and two pinholes before their power can be measured 
independently. The loading procedure is described in section IHl 

Figure [7| shows results for the Rabi spectrum with a 3 ms long tt pulse. A constant 
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Figure 8. Line center {ujc, A) and RMS widtli {a, •) of Rabi spectrum for trapped 
atoms as a function of compensating beam power, for a 3 ms pulse, (ctfl ~ 110 Hz is 
the Fourier limited a). The spectrum width is minimized to a Fourier- limited value, 
at a compensating beam intensity which corresponds also to a minimal shift from Whf- 
The power is normalized to the measured value at which the best compensation is 
achieved, = 25 ± 10 nW. 

background resulting from spontaneous F-changing Raman scattering [10] is subtracted. 
The spectrum of free-falling atoms (figure [7^) shows no inhomogeneous broadening and 
a RMS width, a, which is Fourier limited to 110 Hz. A shift in the peak frequency 
(—756 Hz), and a broadening of the line (to cr = 320 Hz) are seen in the spectrum 
of trapped atoms (figure U}p), in fair agreement with the calculated trap depth and 
atomic temperature. This inhomogeneous broadening is not significantly affected by 
the duration of the pulse. The addition of the weak compensating beam, nearly cancels 
the broadening of the spectrum, as well as its shift from the free-atom line center (figure 

Efc). 

Figure |H1 shows the measured RMS width and shift of the trapped atoms as a 
function of compensating beam power, again for a 3 ms microwave pulse. The spectrum 
width is minimized to the Fourier broadening limit at a compensating beam power which 
corresponds also to a minimal shift from the free-atoms line center. 

Figureinishows the measured spectrum for trapped atoms for a 25 ms long microwave 
pulse. A measurement of free atoms with this pulse length is not possible in our setup 
since the atoms fall due to gravity, and leave the interaction region. A Fourier limited 
a = 13 Hz is measured, representing a 25-fold reduction in the line broadening, as 
compared to the line broadening of trapped atoms. A similar measurement with a 
50 ms pulse shows a nearly Fourier limited width (50 times narrower than the trapped 
atoms spectrum), at the expense of a much larger spontaneous photon scattering and 
hence a smaller signal (All four \F = 2,mp ^ 0) states are populated and contribute to 
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Figure 9. Rabi spectrum of the hyperfine splitting of opticaUy trapped *^Rb, with a 
25 ms pulse. A Fourier limited ct « 13 Hz is measured. 



the spontaneous Raman scattering background, which is hence 5 times larger than that 
of an ideal 2-level system). For even longer measurement times spontaneous photon 
scattering prevents further narrowing of the line. 

We measure the spin relaxation rate ^D] to be ~ 3 x 10^^ s^^ for atoms in the trap. 
The addition of the compensating beam induces an increase of only ~ 20%, as expected. 

It should be noted, that for certain transitions a laser frequency can be chosen such 
as the relative light field perturbations on the measured spectrum cancel [101 EH 1121 • 
Although simpler than our scheme, in the sense that only one laser is needed, ours is 
a more general method which does not require the existence of a "magic wavelength", 
where the light shift of the transition vanishes. 

7. Energy Selection 

In this section we present a method for increasing the coherence time in an optical 
trap, by narrowing the energy distribution of the trapped ensemble using a microwave 
"pre-selection" . First, all the atoms are optically pumped into F = 3, where they 
are distributed among the different mp states. Then, a weak "selection" vr-pulse, with 
frequency utg and duration tg (and hence FWHM spectral width AuJs/2tt ^ 0.8 /tg, see 
section IH), is applied. As a result, a subset of the atoms in ||), composed of those 
atoms having resonance frequency in the vicinity of Ug, is transferred to the lower 
hyperfine state Since the resonance frequency is a monotonic function of energy, 

this selection can be viewed as an energy selection of atoms in the trap, i.e. an energy 
band is selected around the energy Eg. Next, a strong and short laser pulse resonant 
with the 155*1/2,-^ = 3) |5P3/2,-F = 4) transition ejects all the "unselected" atoms 
from the trap without causing any effect on the "selected" atoms. The end-result of 
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the above sequence is an atomic ensemble with a narrower energy spread, and hence a 
longer coherence time. 

The number of selected atoms is given by 

/■oo 

Ns= / g{E)F{E)P{uJoiE))dE (37) 
Jo 

where g{E) in the density of states in the trap, F{E) oc exp{—E/kBT) is a Boltzmann 
factor, and P is the Rabi transition probability, given in equation |22l and expressed 
here as a function of uJo{E) to stress the fact that the resonance frequency depends on 
the energy of the atom in the trap. With our technique, it is possible to select a narrow 
energy band around any central energy, enabling, for example, to maximize the number 
of selected atoms (for a given energy width) by selecting the energy with the highest 
density of thermally populated states. For example, assume a 3D harmonic trap for 
which the density of states obeys g{E) oc E"^. If the trap is populated with an atomic 
sample at a temperature T, and a selection pulse corresponding to an energy Eg and 
(FWHM) width AEg <^ ksT is used, then the number of selected atoms is given by 
Ns{Es) ^ g{Es)F{Es)AEs. For a given AEg, the selection energy which optimizes the 
number of selected atoms is given by i?"^* = 2A;bT. A simple selection of the coldest 
atoms (e.g. by lowering the trapping potential) will result in a dramatically smaller 
number of selected atoms. The ratio between the number of atoms selected around E"^^ 
and near i^^ = is given by 



N{E = E^^') ^ {2kBTfe-' 1 ( ksTV 
N{E = 0) ^ (AEgf ^ 2''\AeJ ■ ^^^^ 

We perform a proof-of-principle experiment of the energy- selection scheme, using 
a red-detuned Gaussian trap with a waist of wq = 50 fim, a power of P = 120 mW 
and a wavelength A = 810 nm (See section IHl). In order to probe the resulting energy 
distribution, we perform Rabi spectroscopy of the remaining atoms using a 20 ms 
microwave pulse. As seen in figure ^1 the spectrum of atoms selected with a short 
100 /is pulse (which actually transfers the whole ensemble to ||)) has a width of ~ 200 Hz 
(FWHM) in agreement with the calculated width for our trap power and detuning and 
the temperature of the atoms. The spectrum of atoms selected with a 20 ms pulse has 
a width of ~ 60 Hz showing a > 3-fold narrowing in the energy distribution. Shown 
in figure ITUl are the results of selection at two different values of Eg, demonstrating the 
ability to maximize the number of selected atoms, as explained above. 

In a similar way to the "echo" spectroscopy of section^] below, any broadening of 
the narrow selected energy-slice as a consequence of different effects, such as photon 
scattering and trap instabilities, can be very instructive on the effects causing the 
broadening. For example, a change in the trapping potential after the selection pulse 
and before the measurement, will result in a line-shape that reflects the trap's local 
density of states (LDOS) HHj. 
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Figure 10. Rabi spectrum of energy-selected atoms. Full dots: Spectrum of trapped 
atoms with no pre-selection. Empty dots: Rabi spectrum of atoms pre-selected with 
a 20 ms pulse at a frequency of Whf/Stt — 156 Hz. Empty triangles: Spectrum of pre- 
selected atoms with a 20 ms pulse at a frequency of louf/'^t^ — 81 Hz. 
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Figure 11. Rabi spectroscopy of atoms in a trap with wavelength A = 805 nm, with a 
20 ms pulse, (a) For a microwave pulse area corresponding to a 7r-pulse the sidebands 
are not visible, indicating that the matrix elements for the sidebands are very small, 
(b) For a stronger pulse, the sidebands emerge. 

8. Microwave Spectroscopy in a Dipole Trap: A Quantum Model 

The classical model of section |3] fails to completely describe the spectroscopic properties 
of the trapped atoms. As an example, figure ^2 shows that for a powerful and 
long enough pulse, the resonance line in a Rabi spectrum of trapped atoms develops 
sidebands, which are, as we shall show next, a clear evidence of the quantization of the 
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Figure 12. The eigenenergies of the trapped atoms consist of two manifolds (belonging 
to II) and It)) separated in energy by -Ehf- The matrix elements for microwave 
transitions depend on the overlap of the vibrational states. The trap shown is a 
Gaussian trap, and gravity gives a small slope. 



motional levels in the trap. 

In reference jTT], a semi-classical, dynamic model was used, in which the center- 
of-mass motion of the atoms is solved, and the phase is integrated over the classical 
trajectory of the atoms. The total signal is obtained by averaging over the different 
trajectories. The validity of this model is not clear, since it assumes the same classical 
trajectory for particles in different internal states. Such a model ignores the separation of 
trajectories with same initial conditions due to the difference in potentials, and predicts 
longer coherence times than the measured ones ^T]. 

We present in this section a quantum-mechanical treatment of microwave 
spectroscopy in an optical dipole trap. The existence of a difference in potential for 
atoms in different hyperfine states, means that the external (center of mass) potential 
depends on the internal (spin) state, hence the internal and external degrees of freedom 
cannot be separated and the entire Hamiltonian including both internal and external 
degrees of freedom has to be considered: 



n = H^\i) {i\ + {H^ + E^,)\]) (Tl 

^2 



p2 



li> Ul 



ITXTI, (39) 



where Vi and are the external potentials for atoms in states |J.) and |t), respectively. 
These potentials include the gravitational potential, equal for both states, and the dipole 
potential, which can be written as and = (l + e)?7j, where e = ojhfI^ can be called 
the "perturbation strength" , typically 10~^ — 10^^ in our experiments. We assume that 
the potential can be approximated as an harmonic one. Since E^^ is much larger than 
both Vi and V|, the eigenenergies of the above Hamiltonian consist of two manifolds 
separated by E^f^ and composed of "vibrational" levels with a separation of hjJosc 
(where Uosc is the oscillation frequency in the trap) jl- We enumerate the eigenstates 

tt We use here one-dimensional notations, since the dynamics in our trap is separable. In section [TSI 
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of the lower manifold using regular numbers and those from the upper manifold using 
primed numbers (see figure IT^ . 

When a microwave field Vmw with frequency close to cUhf is applied, transitions 
between the eigenstates of the Hamiltonian corresponding to different internal states, 
are drivenff. The microwave field acts only on the internal state of the atoms, hence the 
matrix elements for the transitions can be written as the free-space matrix elements for 
the internal state transition, times the overlap between the initial and final vibrational 
eigenstates, which correspond to different Hamiltonians: 

{m, i \Vmw\ n', T) = ("^ I n') ■ (| iKiwI T) = ("^ I n') V^^ 

(m', t |Vmw| n, I) = {m' | n) ■ (t |Vmw| i) = {m' \ n) V^i 

{m', t |Vmw| n, t) = (m' | n') ■ (t |Vmw| T) = 

(m, i I^Mwl n, i) = (m I n) ■ (i \V^^\ j) = 0, (40) 

where we have used the orthogonality of thejln)} and states. 

Now we can understand the result of figure ^2 Non- vanishing matrix elements for 
transitions to vibrational states with n' ^ n show up as sidebands in the microwave 
spectrum, if the pulse is long and weak enough so that the power broadening is smaller 
than the typical level spacing. For a microwave pulse area corresponding to a vr- 
pulse (figure ITTk ) the sidebands are not visible, indicating that the matrix elements 
for the sidebands are very small. Only for a stronger pulse (figure ITTb). the sidebands 
emerge. We verified that the sidebands are stronger for smaller values of the detuning, 
and hence larger perturbation and larger (m | n') for m ^ ra', as expected from the 
analysis. The coupling to different vibrational levels, has a strong analogy to the Franck- 
Condon factors in molecular spectroscopy. Note, that the near harmonicity of our trap 
is responsible for the fact that the sidebands are seen even though many levels are 
thermally populated. 

We start by considering a quantum state written as the external product of an 
internal state and a state which represents the "external" degree of freedom. If the 
atoms are initially prepared in their internal ground state ||), their total wavefunction 
can be written as \1/ = ||) (S> or, shortly, In the presence of the microwave field, 

we need to solve the time-dependent Schrodinger equation given by 

[7^ + K,w]^ = ^^, (41) 

with the Hamiltonian of equation EHl We choose to use ^ = 1 to simplify the equations. 
As is usually done in time-dependent perturbation theory, we choose to work in the 
basis of eigenstates of the Hamiltonian without the microwave field (equation inSl) • A 
general state can be written as 

^ = ^a„|n,i) + ^a„/|n',t). (42) 

n n' 

and references ^1 we treat non-separable (chaotic) dynamics. 

ft Since the size of our trap (~ 50 /im) is much smaller than the microwave wavelength (^ 10 cm), the 
momentum of the microwave photon can be neglected (Lamb-Dicke regime |19|) 
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Introducing this into equation |^ results in 

^ aniEn + Kiw) 1^, i) + X] ^"'(^"' + + Vmw) \n', T) 



\n, I) + ^dn' \n', t) 



(43) 



Two equations are obtained by projecting the above equation onto (m, || and (m', ||: 

a„ (m, I |Vmw| i) + «n' {m, j | Vmw| T) (44) 

n n' 

(Em' + E^p) + ^ a„' (m', t |Kiw| T) + ^ 

n/ n 

Using the matrix elements from equation EHJ these equations can be written as 

n' 

{Em' + -Ehf) + ^ Ctn ("^' I n) V||. (45) 

We assume now a monochromatic, linearly polarized, microwave field. Then we can 
write 

V^^ = [^Mw exp{-iujMwt) + ^Mw exp(iu;Mw^)] , (46) 

where fluw is the microwave field Rabi frequency. It is now useful to define a new 
amplitude, which "rotates" at the field's frequency, bn' = a„' exp(zco'Mw^)- In this new 
frame the equations are 

idm = amEm + ^ {m I n') exp(-ia;Mw^) 

n' 

{Em' + -Ehf - t^Mw) + ^ fln {m I n) V|| exp(it^Mwt)- (47) 

n 

When equation HBl is introduced into equations HH the resulting expression includes 
terms with a exp(— 'i2co'Mw''^) time-dependence. These terms oscillate so fast, compared to 
every other time variation in the equations that they can be assumed to average to zero 
over any realistic time interval, and can be neglected (rotating- wave approximation). 
This results in 

idm, = OimEm 2^l^w ^ ^ ('^ | ) ^n' 

n' 

ibm' = bm'{Em' + Amw) - f^MW ^ ("^' I ^) ^n, (48) 

n 

where Amw = -^hp — "-^mw is the microwave detuning. 
9. Ramsey Spectroscopy in a Dipole Trap 



In general, each of the eigenstates in the ||)-manifold is coupled to many eigenstates 
in the ||)-manifold by the microwave field, hence there is no simple prediction of 
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the wave function after microwave irradiation. However, in order to calculate the 
outcome of a Ramsey experiment (consisting of two microwave pulses of short width 
t, separated by a much longer "dark" period of duration r), we only need to solve 
equations HHl in two limiting cases: first, the short and strong limit (i.e. t <^ uj~J^^ 
and f^Mw ^ Amw + {En — En') for all n,n'), and next, the free evolution limit (where 

In the simple limit of short microwave pulses equations EHl reduce to 

idm = -^J^Mw ("^ I '^O bn' (49) 

n' 

ibm' = -^^Mw ^ {m \n)an. (50) 

n 

Both {\n)} and {|^')} are complete basis that span the vibrational part of the wave 
function. Hence, we can express the vibrational part of \n',^) using {|^)}. The 
coefficients in these expressions, denoted Cn, are given by c„ = J2n' I Using 
also the fact that 6„/ = J2n I ^) ^n, and some algebra, it can be shown that 

These last equations show that there is no coupling between and c„ for m n, 
hence the solution of the coupled equations for every m is given by equation |22l This 
means that if we start with a state = (S> the state after an on-resonance 

pulse of duration t is given by {ip) ® [cos {flMwt / |i) + i sin {fl^^t / 2) If)]. Hence, for 
a TT pulse the vibrational part of the wavefunction is only "projected" on the other 
potential and will result in i ||,^). A ^-pulse will result in the coherent superposition 
state ^(|i,V^)+^ IT, V^)). 

For the "dark" periods we go back to equation HHl using Qmw = 0, hence 

^'^m ~ (^mEm (53) 
ibni' = hm'{Em' + Amw), (54) 

which means that the time evolution of I"?/", t), and \ip-,i) is given by 

T) = exp [-i{H^ + Amw) t] 1^(0), T) (55) 
= exp H/Z^t] 1^(0), i). (56) 

We are now prepared to calculate the result of a Ramsey experiment. Let us 
assume that atoms are initially prepared in their internal ground state, and also that 
their vibrational part is an eigenstate of Vi characterized by the quantum number n. 
Their wavefunction can then be written as {ipn, I), where we use the notation ipn instead 
of \n) to stress that \ipn,^) is not equal to |^, t)- The atom is irradiated with a short 
f -pulse to generate the wave function {\ipn, i) + i li^n, ^)) ■ After some time r this 
state will, in the rotating frame of the microwave field, evolve into 

1 i 

exp [-iHir] |^„, i) + ^ exp [-i {H^ + Amw) t] {tpn, T) • (57) 
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Then the atoms are irradiated with a second |-pulse generating the wave function 
^ {exp [-iHir] - exp [-i {H^ + Amw) t]} {ipn, i) 



2 



+ - {exp [-i {H^ + Amw) r] + exp [-tH^r] } t) . (58) 
The population of state |t) is then, 

which can also be written as 

Pt = i |l + Re eiH,+AMw)r^-iH,r ^^^j | _ ^gQ) 

The expression inside the square brackets is, in general, a complex number which can 
be written as Aexp(2AMwT + 4>n), where 0n(T) is a slowly varying phase depending on 
the initial state \ipn) and the dynamics of it in the trap. Then we can write 

Pt = i (1 + \{^n |e*^i^e-*^T-| cos [Amw^ + <pn (r)]) . (61) 

Equation inU shows that, when Amw is sufficiently large compared to variations 
in 4>n (t), scanning cUmw for a fixed r yields the usual Ramsey fringes with a contrast 
given by \{ilJn |e*^^'^e~*''^T^| ■i/'n)!- The fringes contrast is actually the "fidelity" of the 
external motion in the trap, since it can be viewed as the overlap between a "desired" 
state e~*^J^^ \ipn) and the actual state in a perturbed environment e~^^^'^ iV^n)- The 
fidelity was first proposed by Peres [IHj as an indicator of the stability of a quantum 
system, in an analogue way to characterizing the stability of a classical system with the 
Lyapunov exponent jUj- Alternatively, the Ramsey fringes contrast can be interpreted 
as a Loschmidt echo |1H]; that measures the overlap of a state evolved forward in time 
(with ifj), and then backward in time with a perturbed hamiltonian Hp Finally, since 
the initial wavefunction is an eigenstate of Hi, the contrast of the Ramsey fringes can 
also be written as a time-correlation function {{ip {t = 0) \ {t = r))|. If Vj (x) = (x) 
then an eigenstate of H^ is also an eigenstate of H^, and clearly this contrast equals 
unity. In general, if we start with an eigenstate of H^, the projected state will not be an 
eigenstate of H^, and therefore will evolve in the new potential. Such quantum dynamics 
causes the overlap {ip [t = 0) \ {t = t)) to decay in an interesting way, which depends 
on the type of the underlying classical dynamics (being regular, chaotic or mixed) and 
the strength and type of the perturbation [44]. 

Such quantum dynamics and in particular the decay of fidelity or Loschmidt echo 
in chaotic systems has been the topic of intense theoretical and numerical studies in 
recent years (see for example references |igEniEIlElESlElE3EniEIlEHlEi), 
mainly because fidelity is also the standard measure for loss of information in quantum 
computation jnOj- However, experimental studies of chaotic and mixed systems are 
still a missing chapter, since they require the preparation of highly-excited, pure 
quantum states to avoid "spreading" these interesting effects by averaging over an 
inhomogeneously broadened system. 
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Figure 13. Ramsey spectrum of trapped atoms, measured for t = 1 ms (full line) and 
T = 5 ms (dashed line) . The ensemble-averaged fringe contrast decays rapidly because 
the cosine terms from different populated states get out of phase. 



This difficulty is clearly manifested in our experimental system, composed of a 
thermal ensemble of atoms incoherently populating more than 10^ eigenstates (as 
opposed to an initial single vibrational eigenstate, considered so far). The total 
population in ||) is now given by an average of P| over the initial thermal ensemble. 
If we assume that \{ipn |e'^J-'^e~*^T^| ipn)] = C does not depend on n, then the total 
population is given by 

A = I (^1 + C F{n) cos [AmwT + (pn (r)] j . (62) 

where F{n) = exp(— i^n/fc^T)/ ^^exp(— i^^/Zc^T) is a Boltzmann factor. 

We have shown in reference ^3] that in the small-perturbation regime, achieved 
experimentally for large values of the detuning, an eigenstate of (in the lower 
manifold) is coupled by the microwave field mostly to the corresponding eigenstate 
of (in the upper manifold), and hence the fidelity of each state is nearly unity. 
However, since 0„ (r) = {En — En'=n)t depends on the initial state, the ensemble- 
averaged fringe contrast in a Ramsey experiment will decay rapidly. For this small- 
perturbation regime, the system acts as an inhomogeneously broadened ensemble of 
noninteracting two-level systems and thus the Ramsey-fringe decay time can be simply 
estimated by 1/ {2 Arms), where Arms is the RMS spread of the resonance frequencies 
for |n) — >■ \n' = n) transitions, taken over the thermal ensemble. 

As an example, figure El shows Ramsey fringes of trapped atoms, as measured 
in our experimental setup. The atoms are loaded into a 50 mW horizontal laser beam 
focused to a 1/e^ radius of 50 /xm, and with a wavelength of A = 800 nm, yielding a 
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trap depth of Uo/kfyT = 1.5 (See section |2)). Under these conditions the fidehty of each 
eigenstate as nearly unity. Nevertheless, the contrast of the Ramsey fringes at r = 5 ms 
is reduced to ~ 1/3 of its value at r = 1ms. In Fig^l this measured Ramsey fringe 
contrast is shown as a function of time. As seen, the Ramsey fringe contrast decays on a 
time scale of 2.4 ms, due to the variation of over the thermally populated states. This 
results are in agreement with a calculated decay time of 2.7 ms, assuming a thermal 
ensemble in a harmonic trap clipped at 1.5 ksT. 

In conclusion to this section, Ramsey spectroscopy of optically trapped atoms 
occupying a pure state can yield important information on their quantum dynamics 
and in particular on the fidelity (or equivalently the Loschmidt echo) as they are 
perturbed by the small and well controlled difference in optical potential acting on 
the two internal states. However, in most experiments and in ours in particular, pure 
and highly excited states are practically inaccessible, and for thermal ensembles the 
interesting quantum effects are overwhelmed and completely smeared out by the rapid 
inhomogeneous dephasing of the system. In the next section we show how an echo-like 
scheme can suppress this inhomogeneous dephasing and enable us to directly measure 
the quantum dynamics of thermal ensembles of optically trapped atoms. 

10. Echo Spectroscopy 

The results of the previous section indicate that the decay of the Ramsey fringe contrast 
is not only a fingerprint of decoherence (an irreversible effect) but also a consequence 
of dephasing. Dephasing, which causes the ensemble averaged signal to decay, can 
be reversed, at least partially, by stimulating an effective "time reversal", as has been 
reported for spin echoes and photon echoes jHSESj, and more recently for a motional 
wave packet echo using ultra cold atoms in a one-dimensional optical lattice jHl]. We 
achieve such reversal in dephasing by adding a 7r-pulse, which inverts the populations 
of II) and It), between the two | pulses jUHl- If the vr-pulse is exactly in the middle 
between the two |-pulses the two parts of the superposition state generated by the 
first ^-pulse spend an equal amount of time in both levels, and are therefore exactly in 
phase with each other at the time of the second |-pulse. This means that a "coherence 
echo" appears at the time of the second |-pulse even for a system that has dephased 
completely before the 7r-pulse. The coherence echo is observed by seeing that all the 
atoms return to the initial state after the ^-vr-^-pulse sequence. If the time ti between 
the first | and the vr-pulse is kept constant, and the time T2 between the vr and the 
second |-pulse is swept, the coherence echo is seen as a decrease in the population of |t) 
for T2 = Ti (see figure IT^ . The population for T2 = ti is referred in what follows as 
the "echo signal" where = indicates full revival of coherence and = ^ indicates 
complete dephasing. 

A similar calculation as the one in the previous section shows that the echo signal 
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Figure 14. Echo signal (P| ) measured as a function of the time T2 between the 7r-pulse 
and the second ^-pulse, for a fixed ri. A dip in P| is seen, showing a "coherence echo" 

at T2 = Tl . 

for an initial state {ipn) is now 

A = ^ [1 - Re ((^„ I e'"^^ e^^T- g-^^^^ 6"*^^- | . (63) 

The main difference between tlie echo signal (equation and the Ramsey signal 
(equation is that for the echo P| no longer depends on Amw Hence, when measuring 
a thermal ensemble, the condition {n' \ n) ~ 5nn' for all initially populated vibrational 
states now ensures that P| ~ for all r jUj. In other words, after dephasing for a time 
r a TT-pulse stimulates a coherence echo of the ensemble averaged signal at time 2r. 
Some intuition can be gained, by switching to a new basis {v?a;}, defined by 
= 0)) = exp(— iiZ-fr) !■?/;„). Now, the echo signal can be written using a time 
correlation function: 

P| = ^ {1 - Re [exp (zE^t) {^t = 0) | Mt = r))]} , (64) 

where \ipkit = r)) = exp(— iifj^r) |(y5fc(t = 0)). 

The results of echo spectroscopy with the same trap parameters as used before for 
Ramsey spectroscopy in section IHl are also presented in figure UHl For these parameters 
{n' I n) ~ 6nn' is a good approximation for all thermally populated states. As explained 
before, we subtract from the signal contributions to the population of |t) due to F- 
changing Raman transitions induced by the trap laser and normalize to the signal after 
a short vr-pulse, which transfers the whole population of ||) to If). This corrected signal 
is denoted P|. A coherence echo (P| << 1/2) is clearly seen long after the Ramsey 
fringe contrast has decayed. On a time scale of ~100 ms the echo coherence decays, for 
reasons that will be analyzed further in section El 
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Figure 15. Ramsey fringe contrast (o) and echo signal P| (x) measured as a function 
of the time between the two ^-pulses. For the echo signal the value represents 
complete coherence, and the value 1/2 represents complete dephasing. The trap laser 
wavelength is 800 nm. A coherent echo (P| << 1/2) persists long after the Ramsey 
fringe contrast has decayed. 



11. Quantum Dynamics 

Since echo spectroscopy is seen to cope with the inhomogeneous broadening, we can 
use it to observe the fingerprints of the fidehty even for a thermal ensemble. When 
{n' I n) 7^ 6nn' a good echo signal is no longer expected, since each vibrational 
state is coupled to several vibrational states by the microwave fields, and therefore 
\{^^{t = 0) I = r))| < 1. 

Ignoring gravity, Vj^ and V-^ are just the dipole potentials, and are related by 
^ = (1 + e)Vi, where e = ujjip/6. Hence, the detuning serves as a control parameter for 
the perturbation strength. We perform echo spectroscopy as a function of time between 
pulses for different wavelengths of the trap laser (and hence different perturbation 
strengths, e) while keeping the trap depth constant by adjusting the power of the 
trap beam. The results are shown in figure ^1 For a large detuning (A = 805 nm, 
which means a "weighted detuning" 6 = 16.7 nm and hence a perturbation strength 
e ~ 4 X 10"'') a good echo (P| << 1/2) is seen independent of the time between pulses, 
as also seen in figure IT^ For an intermediate detuning (A = 798.25 nm, e = 9 x 10~^) 
damped oscillations to a level smaller than 1/2 are seen, and for a small detuning 
(A = 796.25 nm, e = 2 x 10~^) a complete decay of the echo coherence (P| = 1/2) is 
observed, followed by partial revivals at later times. 

The interpretation for the large detuning (small perturbation) regime, where 
(n' I n) ~ 6nn', was given above. For an intermediate detuning, a small but 
significant coupling to other {n' ^ n) vibrational levels exists. When atoms are 
transferred from Vj (x) to (x) by the microwave field, a parametric excitation of 
the atomic wavepacket is induced, thereby exciting "breathing" modes. However, the 
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Figure 16. Echo signal Ff measured for three different trap laser wavelengths A. 
Lower curve: A = 805 nm. A good echo signal (P| << 1/2) is seen almost independent 
of time between pulses. Middle curve: A = 798.25 nm. The echo signal oscillates 
around a value smaller than 1/2 with partial revivals at r = Tosc/^ and r = Tqsc, where 
Tosc = 3.6 ms is the measured trap oscillation frequency in the transverse direction. 
Upper curve: A — 796.25 nm. After a short time the echo signal completely disappears 
(P| = 1/2), but partly revives again at r = 3.3 ms, close to Tqsc- 
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Figure 17. Numerical calculations of the echo signal in a 2D harmonic oscillator, with 
gravity, having the same oscillation frequency as the experimental data in figure [TBI A 
similar calculation shows that without the effect of gravity, the revival at r = 1.8 ms 
is complete. 



change in the optical potential also changes the gravitational sag, and hence excites 
"sloshing" modes along the vertical axis. A wavefunction parametrically excited in 
an harmonic oscillator will revive after t = m x \tosc (for = 1,2, ...) yielding 
I l^(pn (t = 0) I (t = m X i Tosc) ) I = 1- ^ sloshing mode will revive after r = m x Tosc 
inni EZI • These revivals are seen in the middle curve of figure El as a partial revival 
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Figure 18. Long time level of echo signal as a function of e, the relative difference 
between potentials, (x): Experimental measurement. Solid line: Calculation of the 
ensemble average of | {n' — n \ n)\ for a 2D harmonic potential with oscillation time 
of 3.6 ms. 



in the echo signal at r = 1.8 ms and a stronger one at r = 3.6 ms, exactly equal to 
the measured trap oscillation frequency in the radial direction. For the upper curve 
in figure El a faster initial decay of the echo coherence is seen, as expected for a 
stronger perturbation, followed by a clear revival at r ~ 3.3 ms. The origin of these 
revivals is verified in numerical calculations of the echo signal for typical single states 
in a 2D harmonic oscillator with gravity, having the same oscillation frequency as our 
experimental data (see figure IT7|) . A similar calculation shows that without the effect of 
gravity, the revival at r = 1.8 ms is complete. The lack of a perfect revival at r = Tosc 
in figure Uni is due to the anharmonicity of the Gaussian trap. Note, that our technique 
is sensitive enough to map the quantum dynamics of a system, due to a perturbation 
approximately three orders of magnitude smaller than ksT. 

For a sufficiently long time r the wave packet oscillations of figure EJ damp due to 
a complete dephasing of the dynamics. At such a long time a simple expression for the 
echo signal can be given. In particular, assuming random phases between all vibrational 
states yields the simple relation 

\{cp^{t = 0) \ if^it = t))\ = \{n' = n \ n)\\ (65) 

Substituting this into equation IHH and averaging over the ensemble yields the expected 
long-time echo signal. We performed this calculation numerically for a 2D harmonic 
trap, in gravity, with our measured oscillation frequency and a thermal ensemble with 
a temperature of 20 /xK clipped at our trap depth of 1.5/cbT. The results are shown in 
figureUHl together with the measured long-time echo signals as a function of perturbation 
strengthf. As seen, the calculation for the harmonic trap and the data points for the 
Gaussian trap show the same qualitative behavior, of improved long time echo when e 

f The ensemble average for the clipped thermal distribution in the harmonic trap was calculated 
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becomes small. 

We end this section with a brief discussion of the validity of the "no mixing" 
condition {n \ n') = 6n,n' in our experimental system. We start by noting that typically 
10^° — 10^^ states of our Gaussian optical trap are thermally populated (this is easily 
estimated by noting that we typically trap ~ 10^ atoms at the standard phase space 
density provided by laser cooling techniques, ~ 10~^). However, since the duration of 
the experiment is often much shorter than the longitudinal oscillation time of atoms 
in our trap, this longitudinal motion is essentially frozen, and we can consider only 
the 10^ — 10^ thermally populated transverse vibration states. Perturbation theory 
indicates that a necessary condition for substantial mixing of states is that the energy 
shift of populated states is comparable to the energy spacing between adjacent states. 
This means that for 10^ thermally populated states the relative energy shift needs to 
be much smaller than 10~^ to ensure the "no mixing" condition. In our experiments 
we typically observe mixing (i.e. a fast decay of the echo coherence) only for much 
larger relative perturbations of ~ 10~^. The solution for this "paradox" is that state 
mixing by a perturbation can be drastically suppressed by strong selection rules, making 
the sufficient condition for mixing ~ 10^ times larger than the necessary one. For our 
Gaussian trap such selection rules are provided by the fact that the transverse motion of 
the atoms is completely separable into two independent one- dimensional motions, with 
only ~ 10^ one-dimensional states thermally populated in each direction. Orthogonality 
thus ensures mixing only between motional states in each direction separately, whose 
nearest-state separation is ~ 10^ times larger then for the general two-dimensional 
states. Such x 1000 enhancement in the stability of the quantum states against 
perturbation in our separable system is only one of many similarly dramatic difference 
of the quantum properties between separable and non-separable dynamical systems. 
Further discussion of such differences exists in section^] and in references |44[ I45j . 

12. "Multiple TT Pulses" Sequence 

We turn now to investigate the limitations on the coherence time achieved by the echo 
scheme. According to the time-independent Hamiltonian of equation EHl the echo scheme 
is expected to provide a complete cancellation of the dephasing for arbitrarily long times, 
provided the condition {n' \ n) ~ 5n,n' is fulfilled for all thermally populated states. 
However, in real physical systems this is not the case [HHI- Two types of mechanisms 
cause the echo coherence to decay for long times even when (n' | n) ^ 5n.n'- The first 
type of mechanisms (which we denote as "dynamical" T2 processes) are time-dependent 
processes leading to a time-dependent resonance frequency of the two-level system. This 

as: E (K^m I "m)| I <^)|)%"'^"'"/''«^/ E e'-^-'/'^'s^ where are eigenfunction of the 

n,m n.m 

harmonic osciUator associated with the horizontal motion and |1), ujj:^ are associated with |2) (the 
spring constant is sUghtly stronger), and are associated with the vertical motion where there 
is both a change in spring constant and a change in zero point due to different gravitational sag. The 
sum is over the ~3x 10^ states for which Enm — {n + m) hujosc < Upot- 
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Figure 19. a: Echo pulse sequence, b: "Multiple- tt" pulse sequence. Shown in dashed 
lines are tt pulse which occur at times where phases are compensated, hence they are 
omitted, with the advantage of less coupling to m^? 7^ states, as explained in the text. 
In our experiment an additional 7r-pulse (not shown in the figure) is used immediately 
following the ^-pulse at the end of the sequence, in order to measure decoherence as 
a deviation from P|=0. 



will cause the two parts of the wave function generated by the first f-pulse to acquire 
different phases during the two "dark" periods of time r, hence causing imperfect 
interference at the time of the second ^-pulse. Examples of such mechanisms are 
fluctuations of the trap depth e.g. due to noise in the trap laser power [UH], fluctuations 
in the bias magnetic field giving rise to a fluctuating second order Zeeman shift, and 
spontaneous Rayleigh scattering of a photon from the trap laser. Rayleigh scattering of 
photons does not lead to instantaneous loss of coherence, as does a Raman scattering 
event JH] • Nevertheless, the recoil energy acquired by the atom can significantly change 
its vibrational level, and therefore its resonance frequency. The result is, again, a time- 
dependent resonance frequency of the two level systems. Other heating mechanisms, 
such as pointing instability of the trap laser beam, typically involve a much smaller 
energy change than Rayleigh photon scattering, hence they induce a much longer time 
scale for dephasing. The second type of mechanisms for decay of echo coherence (Ti 
processes) relates to the finite lifetime of the internal states of the atoms, which is limited 
mainly by transitions induced by the trap laser light. 

As stated above, dynamical T2 processes do not cause instantaneous decoherence, 
but are characterized by a typical time scale over which a substantial phase difference 
evolves. If the time between pulses r is larger than this time scale, we expect the coherent 
signal to disappear. The time r can be reduced by adding more (equally spaced in time) 
TT-pulses between the two | -pulses (see figure HH). If the decay of the echo coherence is 
dominated by dynamical T2 processes, and the time scale for variations in the resonance 
frequency is longer than r, we expect a coherent signal to reappear The method has 
a strong similarity with the "decoupling" method used in quantum information schemes 
in NMR, where a repeated fast refocusing completely stops the dynamics (See reference 
[70] and references therein). 
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Figure 20. • : Coherence signal (P|) for tlie eclio pulse sequence (figure EJi) as a 
function of TTotai . U : P-^ for a pulse sequence with six 7r-pulses (figure I19b ^ as a 
function of TTotai- It is seen that adding more 7r-pulses increases the coherence time, 
but also leads to an initial small and rapid partial coherence decay. In the inset the 
short time signal for a pulse sequence with ten 7r-pulses is shown. Wave-packet revivals 
are seen for TTotai = 7 and 14 ms. A : P| for a tt-tt pulse sequence as a function of 
TTotai- A monotonic increase in the signal is seen, due to transitions between m^-states 
within the same F-manifold. 



In this experiment the dipole trap consists of a 370 mW hnearly polarized horizontal 
Gaussian laser beam focused to a 1/e^ radius of 50 /im and with a wavelength of 
A = 810 nm, yielding a trap depth of Uq ~ 100 /xK. The transverse oscillation time 
of atoms in the trap is measured by parametric excitation spectroscopy to be 1.4 ms [7T] 
ensuring (n' | n) ~ 6n,n', for all thermally populated transverse vibrational states. The 
bias magnetic field is stronger (250 mG) and turned on after the atoms are loaded into 
the trap. 

As seen in figureEHl the echo signal starts from P| = (indicating perfect coherence) 
and monotonically approaches -Pt ~ i (indicating complete loss of coherence). The 
(1/e) coherence time Tc is seen to be Tc = 26 ms. This coherence time is already ~ 5 
times longer than the Ramsey fringes decay time, showing the strong suppression of 
dephasing provided by the echo technique \. Next, we add more vr-pulses using the 
pulse sequence shown in figure E| First a |-pulse creates a coherent superposition 
state of II) and |t). After time r the first vr-pulse is applied, followed by the rest of 
the vr-pulses at time intervals 2r. A second f-pulse is applied at time r after the last 
TT-pulse, and followed immediately by an additional vr-pulse (not shown in the figure) 
in order to have P| = for a coherent signal for the even number of vr-pulses that we 
use. The signal as a function of the total time between the first and the last pulse, 
TTotai , for a pulse sequence containing six vr-pulses is also shown in figure EDI As seen 

I The trap in the experiment reported here is ^ 4.5 times deeper than the trap used in section [TUl 
hence the echo coherence time is shorter. The Ramsey decoherence time, which is determined mainly 
by the temperature of the atoms, is almost equal. 
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Figure 21. (•) Coherence time as a function of tlie number of 7r-pulses. A maximum 
coherence time of ~ 65 ms is observed, as a consequence of a trade-off between 
suppressing r2-processes and the increased mixing to different vibrational states. 
(a) Twice the slope of the intermediate regime, corrected for contributions from F- 
conserving transitions. 



in the graph, at long times P| exceeds 1/2, which corresponds to total decoherence for 
a purely 2-level system. We attribute this to a weak population mixing between mis- 
states within the same F-manifold as a consequence of spontaneous Raman scattering 
and off-resonant stimulated Raman transitions involving two photons from the trap 
laser. The net effect of these transitions is to increase the population of the F = 3 
manifold after the subsequent microwave pulses. To isolate and directly measure this 
population mixing effect, we measure the population of F = 3 as a function of time 
between two population-inverting vr-pulses §. From the results, also presented in figure 
EDI we measure the F-conserving transition rate to be 1.2 s^^. Combining these results 
with the (independently measured) 0.6 s~^ rate of F-changing transitions piH l29j. yields 
a population decay rate of 1.8 s~^ for the mF=0 states. 

Since the value of at long times is greater than 1/2 we conservatively define the 
coherence time as the time TTotai between the two f -pulses where P| reaches a value of 
P| = |(1 — i). It is seen that Tc = 65 ms, clearly showing that the additional vr-pulses 
substantially increase the coherence time. 

As shown in figure EDI the multiple- vr sequence suffers from an initial small and 
rapid partial loss of coherence. This is because {n' \ n) is not strictly a delta function, 
hence when more vr-pulses are added, the mixing to other vibrational levels is increased 
and the dynamic effects discussed in section [TT] appear (The asymptotic value of due 
to this mixing effect is |[1 — {n' \ ra)^^"""^^-*] where is the number of vr-pulses. The 

§ To further support our interpretation, we verify that a reduction in the transition rate for larger 
values of the bias magnetic field is observed. Moreover, we verify that the asymptotic slope of the 
multiple- TT echo signal is indeed half the slope of the tt-tt measurement, as expected from the fact that 
the mp population difference of the former is roughly half that of the latter. 
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effect is too weak to be observed for = 1). Shown in the inset of figure 1^ is the short 
time coherence signal for a sequence with 10 7r-pulses. Wave packet revivals appear for 
TTotai =7 and 14 ms (i.e. a time between vr-pulses of 0.7 and 1.4 ms) as expected for an 
harmonic trap with our measured transverse oscillation time of 1.4 ms. 

In figure 1^ the coherence time Tc is shown as a function of the number of vr-pulses. 
It is seen that the coherence time initially shows linear dependence on the number of 
TT-pulses and then reaches a maximum value of ~ 65 ms, a 2.5 fold improvement as 
compared to the simple echo coherence time. The maximum coherence time is given 
by a trade-off between suppressing the dephasing due to T2 processes, by adding more 
TT-pulses, and the increased dephasing due to mixing to other vibrational states, that 
the additional vr-pulses induce. 

The decay rate after the initial fast decay shows the residual decoherence, not 
suppressed by the multiple- vr scheme. We fit the slope of the intermediate, nearly linear 
section of the coherence decay curves with a straight line, and correct for the long time 
slope (We show in figure EH twice the slope since the decay of coherence contributes 
only 1/2 to P|). We see that the addition of vr-pulses improves the intermediate slope 
by a factor of 6. However, as explained before, the initial mixing to other vibrational 
levels prevents us from fully exploiting this improvement. 

13. Quantum Dynamics and Chaos 

The coherence properties of trapped atoms are related to their dynamics in the trapping 
potential. Loosely speaking, the fact that the same electron which is probed by the 
spectroscopy is used for the trapping, means that there is an inherent interplay between 
dynamics and spectroscopy. Already at the first demonstration of hyperfine spectroscopy 
in a dark optical trap [TT], it was conjectured that dynamics and spectroscopy are 
coupled in dipole traps, and that the classical trajectories of the atoms in the trap affect 
the coherence time of a spectroscopic measurement performed on them. This interplay 
between dynamics and spectroscopy means that precision spectroscopy in optical traps 
can serve also as a sensitive probe of the dynamics. 

The quantum manifestations of different types of classical dynamics is still an issue 
of unsettled debates, since the unitary evolution characteristic of quantum mechanics, is 
inconsistent with the classical definition of chaos, i.e. an exponential sensitivity to the 
initial conditions. A widely accepted signature of quantum chaos is the correspondence 
between the spectral statistics of quantized classically chaotic systems and the one of the 
canonical ensembles of random matrix theory [72] . An alternative approach is to look at 
the stability of the dynamics with respect to small changes in the Hamiltonian, as first 
proposed by Peres jlH] . In this framework, the "fidelity" denotes the overlap between a 
state evolved by a Hamiltonian with the same state evolved by a slightly perturbed 
Hamiltonian Hi [ISj. Despite the renewed theoretical and numerical interest in the 
fidelity and its properties (see references ggEUlEIlElllSSlElESlEniEZlEiEni), 
experimental studies in chaotic and mixed systems are still lacking, mostly due to 
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the difficulty of preparing highly-excited pure quantum states. The use of "echo 
spectroscopy" eliminates the need for pure quantum state preparation, since quantum 
fidelities can be experimentally observed even in a thermal ensemble and for extremely 
highly excited quantum states. 

As a specific example, we use in reference p3| "echo spectroscopy" to experimentally 
measure quantum dynamics of ultra cold *^Rb atoms trapped in atom-optics billiards 
|731 17^ 17^ I76j ■ with underlying chaotic or mixed classical dynamics. In this experiment 
the trap is a light-sheet wedge billiard, made from two crossed blue detuned light sheets 
defining the billiard walls and where gravity confines the atoms in the vertical direction 

The light sheets have in-focus (1/e^) dimensions of 20 x 250 fim, and by the use of 
cylindrical lenses mounted on rotational stages, the wedge angle are adjusted in order 
to control the classical dynamics||. The confinement in the longitudinal direction is 
provided by the diffraction of the beams. The very elongated shape of the trap allows 
us to consider only the transverse motion and neglect the longitudinal one, which has a 
timescale much longer than the experiment time. The temperature of the atoms is much 
larger than the mean level spacing in the trap, and the atoms typically occupy many 
(up to ~ 10^) states in the trap. The measured echo signal is the ensemble average of 
all of them. 

As shown in reference 03], the echo amplitude can be expressed as a function of 
the local density of states (LDOS). The LDOS denotes the local average of the absolute- 
value-squared matrix elements of the transformation matrix from eigenstates of to 
eigenstates of H^. Formally, it is simply | (m^ \ni)f as a function of — averaged over 
a ensemble of neighboring with approximately the same energy. The width over which 
the LDOS is nonvanishing is denoted the "bandwidth", and if it is large, we expect a 
rapid decay of the echo coherence. In the perturbative regime the LDOS displays system 
specific features, despite the fact that the underlying classical dynamics is chaotic. In 
particular, it is evident from the semiclassical calculations of E] that for atom 
optics billiards, where the inherent perturbation is localized on the billiard walls, the 
LDOS has pronounced peaks for — corresponding to En — Em = h/ru, where tu is 
the typical time between encounters with the wall (~ 15 ms in our experiments). This 
means that the echo signal will show partial revivals for r = tu |44j . 

In figure |221 the decay of the echo signal for different perturbations is presented 
for a wedge with vertex half-angle a = 52.5° (see inset) where the dynamics is fully 
chaotic. For small perturbations a nonmonotonic decay is seen with a partial revival of 
correlations for r ~ r^i as predicted in |3^. The revivals are seen despite the fact that 
due to the "high" temperature a large band of energies is occupied in the trap. Partial 
revivals of correlations in traps where the dynamics is separable are generally expected at 
time scales equivalent to the ID level spacing, but their surprising observation here, for a 

II The structure of phase space in an ideal wedge biUiard can be tuned from stabihty to chaos by 
varying the vertex half-angle, a [77j. For a < 45° phase space is mixed, and the size of the stable 
islands oscillates as a function of a, a behavior dubbed "breathing chaos" . For a > 45° the system is 
fully chaotic. 
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Figure 22. Echo signal for a liglit slieet wedge witli cliaotic classical dynamics, 
for different perturbation strengths (trap laser wavelength tuned from A = 775.9 to 
A = 779.7). ♦: e = 1.44 x lO^^, V: e = 1.90 x IQ-^, A: e = 2.43 x lO'^, o: 
e = 3.80 x 10~^ and ■: e = 1.52 x 10~^. For small perturbations a nonmonotonic 
decay with revivals around r — tu is seen, whereas for large ones a monotonic decay 
is observed. Inset: CCD image of the trap laser at the focal plane, showing a wedge 
with a = 52.5°). 



trap with chaotic motion, demonstrates that they are a far more widespread phenomena. 

For larger perturbations a crossover to a regime where the decay is monotonic 
is observed. In this regime the decay is independent of perturbation strength. The 
perturbation is large enough so the overlap of equivalent eigenstates is small, and this 
indicates that the effects of quantization of the trap levels should not play a role and 
a classical description might be possible. Since the echo amplitude can be viewed as a 
propagator |44j, and in our system the thermal de-Broglie wavelength is much smaller 
than the billiard's dimensions, it is possible to use a semiclassical propagator f72|. It 
is then seen that the classical trajectories contributing to the ensemble average of the 
echo amplitude, are those that after evolving forward in time in and ifj, and then 
backwards in time in and if^, return to the vicinity of their initial position. These 
trajectories we divide into two types: those that during this propagation hit the wall, 
and those that do not. Since and ifj are highly different mainly in the vicinity of the 
wall, then the action integral along the first type of trajectories will yield a very large 
phase and the contribution from these trajectories to the ensemble average of the echo 
amplitude will average out. However since the second type of trajectories does not feel 
the difference between the potentials it will retrace its forward propagation backwards 
in time causing the action integral to vanish. These trajectories will give a perfect 
contribution. Therefore the echo amplitude in this regime simply measures the classical 
probability that the particles have not yet hit the wall. Calculations of the echo signal 
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using this model yields good agreement with the experimental measurements. 
14. Summary 

In this tutorial we present our experiments on the microwave spectroscopy of the ground 
state hyperfine splitting of optically trapped atoms. We introduce a number of schemes 
to overcome the inherent perturbation of the atomic levels by the trapping potential. 

First, we present a simple, classical model and use it to analyze the light- 
induced perturbations on optically trapped atoms, and in particular the inhomogeneous 
broadening of the hyperfine splitting, for a number of trap geometries. An "optimal" 
trap, achieved by minimizing the surface to volume ratio of the trap (to almost that of 
a sphere) and by choosing the minimal wall thickness which is allowed by diffraction, is 
shown to allow an atomic coherence time which is ~ 200 longer than existing single-beam 
dark traps, and ~ 1800 times longer than red-detuned traps. 

Next, we demonstrate a scheme for eliminating the trap-induced inhomogeneous 
broadening of the transition, by adding a weak "compensating" laser, spatially mode- 
matched with the trapping laser and with a proper detuning and intensity. Despite 
being tuned close to resonance, this laser induces a negligible change in the dipole 
potential, and does not considerably increase the spontaneous scattering rate. With the 
suppression of inhomogeneous broadening, the atomic coherence time is now limited by 
the much smaller spontaneous scattering time. Whereas in a conventional optical trap 
the ac Stark shift of the line center strongly depends on the temperature of the atoms, 
which may drift considerably, in the compensated trap the suppression of the line shift 
is equally effective for all temperatures. Hence, it provides a means of achieving a higher 
stability of the line center than that achieved by simply stabilizing the trapping laser 
detuning and intensity. For relative spectroscopic measurements, such as the proposed 
measurement of the electron's permanent electric dipole moment (EDM) [THj, only 
stability (and not absolute accuracy) of the line center is of importance. For example, for 
a 10 fiK deep YAG-laser trap, and a compensating beam with a 15 KHz (time averaged) 
frequency stability, locking the relative intensity between both beams to a 1 : 10~^ 
stability, will result in ~ 10-^"^ stability of the microwave line center. 

Conversely, we present a method in which a long microwave pulse is used to select 
a narrow energy band from the atomic ensemble, around any central energy. The rest 
of the atoms are expelled from the trap using an on-resonant laser, without perturbing 
the selected atoms. With our method, the central energy can be chosen to maximize the 
number of selected atoms by selecting the energy with the highest density of populated 
states. 

We also demonstrate that a macroscopic coherence, lost as a consequence of 
dephasing of the different populated vibrational levels, can be efficiently revived by 
stimulating an echo if the trap detuning is large enough so \{n' = n \ n)\ ^ 1 for 
all thermally populated states. This suppression of dephasing due to perturbations 
induced by the trap yields a dramatic increase in the coherence time for trapped 
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atoms, that may find important applications for precision spectroscopy and quantum 
information processing. We furtfier demonstrate tliat ecfio spectroscopy can be used 
to experimentally map the quantum dynamics of trapped atoms, showing a crossover 
between quantum and classical regimes. In this way it can serve as an extremely sensitive 
experimental tool for investigating quantum chaos in atom-optics billiards even when 
more than 10^ states are thermally populated, thereby avoiding the need of preparing 
highly excited pure quantum states. 

We demonstrate that the dephasing in microwave spectroscopy of optical trapped 
atoms due dynamical changes in the trap parameters, or other similar processes, can 
be suppressed beyond the suppression offered by the simple echo, by using an improved 
pulse sequence containing additional vr-pulses. The achieved coherence time is limited 
by increased mixing between transverse states and, to a lesser extent, the lifetime of 
the internal states of the atoms. Both these factors are expected to be substantially 
smaller for a trap laser with much larger detuning, such as in [HH]- The demonstrated 
pulse sequence may also find use in precision spectroscopy of a periodic effect, where 
the TT-pulses can be synchronized with the period of that effect. 

Finally, we show that "echo spectroscopy" can serve as an important tool for the 
study of the properties of the fidelity (or Loschmidt Echo) and its relation to the 
underlying classical dynamics of a quantum system. The main difficulty in experimental 
studies of these matters, is the preparation of highly-excited pure quantum states. This 
difficulty can be tackled by using echo spectroscopy, which eliminates the need for pure 
quantum state preparation and allows the experimental observation of quantum fidelities 
even in a thermal ensemble and for extremely highly excited quantum states. 

A first work in this direction in our experimental system was done in references 
|441 145j . Using a perturbative treatment we show that the coherence of the echo scheme 
is a function of the survival probability or fidelity of eigenstates of the atoms in the trap 
. The echo coherence and the survival probability display "system specific" features, 
even when the underlying classical dynamics is chaotic. We performed echo spectroscopy 
in chaotic and mixed atom-optics billiards, as a function of the perturbation strength, 
and observed two different regimes |45^. First, a perturbative regime in which the decay 
of echo coherence is non-monotonic and partial revivals of coherence are observed, 
as opposed to the prediction from random matrix theory. These revivals are more 
pronounced in traps with mixed dynamics as compared to traps where the dynamics 
is fully chaotic. Next, for stronger perturbations, the decay becomes monotonic and 
independent of the strength of the perturbation. In this regime no clear distinction can 
be made between chaotic traps and traps with mixed dynamics. 

Some of the ideas presented here, have been found useful by other experimental 
groups. A "compensating" technique has been used (7^1 to compensate the Stark shift 
in an ion trap quantum processor. Recently, two other groups working on quantum 
information processing with single ions and atoms have reported the use of coherence 
echoes jSHlinH!. 

The combined possibilities offered by our experimental system (control of dynamical 
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phase-space and spectroscopy) can have impact on other fields of study: First, the effect 
of interactions between the particles trapped in a billiard is an interesting question in 
the context of semiconductor quantum dots. In atom-optic billiards, the interactions 
can be tuned in a controlled way, by changing the density of the trapped atoms or by 
the use of a magnetic field near a Feschbach resonance to change the atom's scattering 
length [HT]. Next, open systems, which are characterized by a scattering matrix (and 
not by eigenstates) can be investigated, and related to results in quantum dots, such 
as universal conductance fiuctuations jH2]- Finally, chaotic atom-optic billiards, which 
their unique ability to move from the quantum to the classical regime, can be used for 
experimental studies of parametrically-dependent Hamiltonians jHB] and driven quantum 
systems jH3j. 
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